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Abstract 
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To better understand the contribution of cosmic ray muons to the CUORICINO 
background, ten plastic scintillator detectors were installed at the CUORICINO site 
and operated during the final 3 months of the experiment. From these measurements, 
an upper limit of 0.0021 counts/ (keV-kg-yr) (95% C.L.) was obtained on the cosmic 
ray induced background in the neutrino less double beta decay region of interest. The 
measurements were also compared to Geant4 simulations. 
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1 1. Introduction 

2 Understanding the nature of neutrino mass is one of the key topics at the frontier 

3 of fundamental physics. One of the best opportunities for investigating this problem 

4 is searching for neutrinoless double beta decay (Ov(3(3), a transition in which a nucleus 

5 (A, Z) decays into a daughter (A, Z + 2) with the emission of two electrons but no 
e (anti-)neutrinos. 

7 The CUORICINO experiment was a 130 Te-based search for 0^/3/3. It consisted of 

a an array of 62 tellurium dioxide (Te02) bolometers with a total mass of 40.7 kg. It 

9 was operated at the Laboratori Nazionali del Gran Sasso (LNGS) in Assergi, Italy, 

10 from early 2003 to June 2008. The CUORICINO detector was built as a prototype 

11 for the CUORE experiment, which will have 19 CUORICINO-like towers and is 

12 presently under construction at LNGS. 

13 The CUORICINO crystals were arranged in a tower made of 13 levels, 11 with 

14 four 5x5x5 cm 3 crystals and 2 with nine 3x3x6 cm 3 crystals. Each crystal was 

15 operated as a bolometer able to detect an energy deposition by recording the resulting 

16 temperature increase with a neutron transmutation doped Ge thermistor [T] . In the 

17 case of Ou/3/3, the summed energies of the electrons and recoiling nucleus would result 
is in a mono-energetic peak at the 0^/3/3 transition energy of 2527.518 ± 0.013 keV for 

19 130 Te 0. 

20 The detector operated at ~10 mK, cooled by a dilution refrigerator and sur- 

21 rounded by several layers of shielding. Directly above the detector was a 10 cm thick 

22 layer of low-activity "Roman" lead (from ancient Roman shipwrecks). Around the 

23 sides of the detector were several layers of thermal shields and a 1.2 cm thick cylindri- 

24 cal Roman lead shield. The thermal shields were made from electrolytic copper and 

25 totaled at least 1.5 cm in thickness. Outside the cryostat was a 10 cm low-activity 

26 lead shield and a 10 cm standard lead shield. The cryostat and shields were sur- 

27 rounded by a Plexiglas box flushed with clean N 2 from a liquid nitrogen evaporator 

28 to avoid radon, followed by a 10 cm borated polyethylene neutron shield. A top lead 
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29 shield was located about 50 cm above the top plate of the cryostat. The entire setup 

30 was enclosed in a Faraday cage to reduce electromagnetic interference. The assembly 

31 is shown in Figure [T] A more detailed description of the detector can be found in 

32 Ref. pQ. 




Figure 1: The layout of CUORICINO showing the tower, the various heat shields, and the external 
shielding 

33 In CUORICINO, any single-bolometer energy deposition in the Ovflfl energy re- 

34 gion is a potential background that can decrease the sensitivity of the experiment. 

35 Cosmic rays are one source of background. The 3200 mwe overburden at Gran Sasso 

36 eliminates the soft cosmic ray component and reduces the flux of penetrating muons 

37 by six orders of magnitude to ~1.1 ^/(h-m 2 ) |3J, with a mean energy of ~270 GeV 

38 and an average zenith angle (8) ~ 35 degrees. The azimuthal distribution reflects 

39 the mountain profile [U [5] . 

40 A muon could produce a bolometer signal by interacting directly in the detector. 

41 Additionally, muons interacting in the detector, shieldings, or surrounding materials 
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42 could create secondary products that might mimic a Ou/3/3 decay. For example, 

43 neutrons produced by cosmic rays are very energetic and thus difficult to block 

44 with shields. Photons emitted in (n, n! 7) or (n, 7) reactions could appear near the 

45 0uf3f3 energy. Neutron production increases with the atomic weight of the material; 

46 therefore, lead shields can be a strong source of muon-produced neutrons. However, 

47 neutron production is mostly associated with showers, so this background may be 

48 effectively identified by coincident events in different bolometers. 

49 Several Monte Carlo simulations have been carried out on cosmic ray-induced 

50 backgrounds but few direct measurements have been made [6l [H [8j [9j [TUl CHI H2J [El 

51 [H]. For the present study, an external muon detector was installed to tag muon- 

52 induced background events in CUORICINO during its last three months of operation. 

53 Section [2] and Section [3] give details of the muon detector setup and performances. 

54 Section [4] is a summary of the Monte Carlo simulations, while Section [5] describes 

55 the data analysis and results. 

56 2. Muon Detector Setup 

57 An array of ten large plastic scintillators placed outside of the Faraday cage, 

58 which surrounds the detector, was used to tag muons. The scintillation counters 

59 were obtained from previous experiments; the various types are described in Table [Tj 

60 The total sensitive surface area of the scintillators was about 3.67 m 2 . A photograph 

61 of four of the scintillators is shown in Figure |2j 



Scintillator 


Length 


Width 


Thickness 


Number 


Label 


(cm) 


(cm) 


(cm) 


of PMTs 


Al 


100 


50 


5 


1 


A2 


100 


50 


5 


1 


Bl 


120 


60 


15 


2 


B2 


120 


60 


15 


2 


CI 


96 


42.5 


3.2 


1 


C2 


55 


64 


3.2 


1 


Dl 


200 


20 


3 


1 


D2 


200 


20 


3 


1 


D3 


200 


20 


3 


1 


D4 


200 


20 


3 


1 



Table 1: Dimensions of the plastic scintillators used 
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Figure 2: Four of the ten scintillators used (types B and D), shown from the top of the CUORICINO 
Faraday cage 

62 The scintillators were deployed to tag as many as possible of the muons hitting 

63 the lead shields while accounting for both the angular distribution of the incoming 

64 muons and the geometric constraints from existing structures. A simple Monte Carlo 

65 simulation reproducing the muon flux measured by MACRO [5] was used to optimize 
ee the placement of the counters. The arrangement of the scintillators is shown in Figure 

67 El 

es Each scintillator was read out by one photomultiplier tube (PMT) attached to 

69 one of its smallest faces, except the type B scintillators which had two PMTs on the 

70 same face with their outputs summed. 

71 The type A and B scintillators were the thickest and were operated alone. For 

72 these scintillators, the energy released by a through-going minimum ionizing particle 

73 was greater than 8 MeV, which was well above any naturally occurring gamma or 

74 beta background as well as most naturally occurring alpha lines; therefore, muons 

75 may be discriminated from background by simply applying cuts on the energy. The 

76 type C and D scintillators were about 3 cm thick and were operated in pairs. For 

77 each pair, one scintillator was stacked on top of the other and a trigger signal was 

78 generated only when they were hit in coincidence (within 120 ns of each other), as 

79 indicated in Figure |4} A 5 cm thick layer of lead was placed between each pair of 
so type D scintillators to further reduce backgrounds. 
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D1-D3 



B1 - B2 D2-D4 




Figure 3: The drawings show the positions of the scintillators around CUORICINO (left: side view, 
right: top view): the large, dark grey box is the neutron shield placed around the detector, and the 
smaller objects arc the scintillators and the top part of the lead shield. The support structures for 
the scintillators have been omitted. 

The signals from the PMTs were sent to the electronics and data acquisition 
(DAQ) systems. The analog electronics stage, constructed from commercial NIM 
modules, was responsible for generating the trigger signals; the analog signals were 
afterwards digitized by a dedicated VME data acquisition system synchronized with 
the CUORICINO DAQ (Figure Q- Each PMT signal was split in two copies: one 
was sent to a threshold discriminator; the other, after being delayed, was fed into 
a VME QDC board (Caen V792 N). The QDC board recorded the charge from the 
PMT (integrated over 120 ns), which was proportional to the energy released in the 
scintillator. The logic signals from the threshold discriminators were also split: one 
copy went to the NIM boards implementing the trigger logic, while the other went 
to a VME TDC board (Caen V775 N). The TDC board recorded the relative time 
between all PMT signals and the trigger, with a nominal precision of 70 ps. However, 
since the typical time resolution of the PMTs was 1-2 ns, the relative time between 
multiple PMT hits associated with a single trigger is known to a few ns, while the 
absolute trigger time is known only to the precision of the CUORICINO DAQ (8 
ms). 
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Figure 4: Principle of operation of the electronics and the DAQ system 



97 3. Detector Operation and Performance 

98 The muon tagging system was operated with CUORICINO from 12 March to 

99 26 May 2008. The system was running ~53% of the time because of CUORICINO 

100 calibrations and downtime for repairs and maintenance; the total live time was 38.6 

101 days. 

102 Figure [5] shows the energy spectrum acquired by one of the type A scintillators. 

103 Two regions are evident: a low energy background region and a broad peak at higher 

104 energies. The low energy background is due to radioactivity, dark noise, and muons 

105 that clip the scintillator, whereas the higher energy peak is mostly due to cosmic ray 
loe muons. 

107 The efficiencies of the detectors were measured above ground in the assembly 

los hall of LNGS, where the muon rate was much larger. The measurement was made 

109 by placing a pair of scintillators (A and B) above and below the scintillator whose 

no efficiency was being measured (C), such that any muon passing through both A and 
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700 



QDC value 



Figure 5: Energy spectrum acquired underground at LNGS using detector Al. The darker region 
(not used in the analysis) is dominated by low energy background, while lighter region primarily 
contains muon events. The X axis is QDC counts, proportional to energy. 
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in B must also pass through C. If Nab is the number of hits occurring in coincidence 

112 in detectors A and B, and Nabc is the number of hits in coincidence between all 3 

113 detectors then the efficiency of detector C is simply rjc = Nabc /Nab- Of course, the 
H4 efficiency depends on the thresholds set by the threshold discriminators; the thresh- 
H5 olds were chosen based on these measurements to be as high as possible while still 
lie maintaining an efficiency close to unity. The individual efficiencies of the detectors 
H7 measured in this fashion were generally greater than 95%; however, there was an 
us additional loss of efficiency from cuts applied in the analysis to reduce background 
n9 (described below). 

120 Muons were discriminated from the background in the thick scintillators (types 

121 A and B) with an energy cut: for these detectors a software threshold was set at 

122 the dip between the signal and background regions in the spectrum (Figure [5]). A 

123 rough estimate of the loss of efficiency due to this cut was obtained by assuming 

124 a Gaussian shape for the muon peak. For the type A (5 cm thick) detector shown 

125 in Figure [5j this cut rejected ~10% of the muons, although for the thicker type B 

126 (15 cm) detectors, the estimated loss of efficiency was less than 1%. 

127 For the thin scintillators, coincidences between different detectors were used to 

128 generate triggers as described in Section [2] and no further cut on the energy of the 

129 events was applied in the analysis, since the muon peak was not well separated from 
no the background in the energy spectrum. 

131 In order to determine the overall efficiency of the setup for tagging muons asso- 

132 ciated with CUORICINO bolometer events, the intrinsic detector efficiency and the 

133 efficiency of the software cuts were combined for each detector. This information was 

134 then included in a Monte Carlo simulation precisely reproducing the CUORICINO 

135 geometry, the positions of the scintillators, and the distribution of the muon flux. 

136 The simulation is described in more detail in Section HI 

137 The total trigger rate of the muon detectors combined was ~14 mHz with no cuts 

138 applied, or ~4 mHz with energy threshold cuts, while the expected signal rate from 

139 the simulation was 1.75 mHz. The difference between the predicted and measured 
uo rate is due to the fact that the trigger thresholds were kept low in order to maximize 
ui the efficiency of muon detection; this resulted in the inclusion of some triggers caused 
U2 by radioactive decays and dark noise. These spurious muon triggers contributed some 

143 background to our measurement through the increased rate of accidental coincidences 

144 between the scintillators and bolometers, which was taken into account in the analysis 

145 described in Section HI 
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146 4. Simulation 

147 Geant4 version 9.^] [15] was used to simulate the muon-induced backgrounds in 
us CUORICINO. The LBE (Low Background Experiment) physics list was used. The 

149 Geant4 capability of event-by-event simulation was employed to follow the whole 

150 sequence of secondary tracks from the initial interaction to the detector, including 

151 the contribution of neutrons generated from muon interactions in the shields. The 

152 complete structure of the scintillators, external shields, internal shields, and detector 

153 geometry was implemented according to the model shown in Figure [3] The propaga- 

154 tion of particles through the rock overburden was not simulated, but was accounted 

155 for as described below. 

156 An external code simulated the muon energy and angular distribution in the 

157 underground laboratory of LNGS. Muons were generated on a 6 m hemisphere in 

158 the underground laboratory according to the angular distribution measured by the 

159 MACRO experiment [5]. The generated muons were then assigned an energy based 

160 on the ground-level energy spectrum for that angle, which was approximated as |16j : 

dN 0.14 • Eq 2 l 7 ( 1 0.054 

dE G i ■ dVt s cm 2 sr GeV V 1 + a E GL cos 6 + 1 + E GL cos 6 

lei where E G l is the energy at ground level, a = 1.1/115 GeV, and (3 = 1.1/850 GeV. 

162 The ground-level energy was then translated into an underground energy based on 

163 the formula [IB] : 

Eu = (E GL + e) ■ e~ bX ^ - e, (2) 

164 where E v is the energy underground, b = 0.4 x 10 -5 cm 2 /g, e = 540.0 GeV, and 

165 X(6, 0) is the thickness times density of the overburden in the given direction. The 
lee advantage of this method is that it includes the correlation between the direction and 
167 energy of the muons underground. The range of above-ground energies simulated was 
lea chosen for each direction such that the underground energies spanned from 1 GeV 

169 to 2 TeV, which corresponds to ~99% of the underground muon flux. 

170 The output of the simulation contained the event number, detector number (scin- 

171 tillator number or bolometer number), hit time, and energy released in the detector. 

172 This output was used to produce spectra and scatter plots, taking into account the 

173 detector response and analysis cuts in order to reproduce the experimental condi- 

174 tions. A Gaussian smearing of 8 keV (full-width at half maximum) modeled the 

175 bolometer resolution. 



*A known bug affecting the neutron inelastic interactions has been fixed in Geant4 9.2: 
: / / geant 4 . cern . ch / support /ReleaseNotes4 . 9 . 2 . html 
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176 The Monte Carlo simulation produced the equivalent of about 3.5 years of data 

177 (~ 8 x 10 6 primary muons). In addition to statistics, the simulations were subject 
us to systematic uncertainties: uncertainty in the primary muon flux and spectrum 

179 (8%) [3J, Geant4 electromagnetic tracking (5%), uncertainty in the muon- induced 

180 neutron yield (40%), and neutron propagation and interaction (20%) [T7]. Analysis 

181 of simulation results will be discussed in Section 15.11 



182 5. Data Analysis 

183 The analysis involved searching for correlations between muon triggers and events 

184 in the CUORICINO bolometer array. A coincidence was defined as a muon detector 

185 event occurring within ±50 ms of a bolometer event. This large window, chosen 

186 based on the time resolution of the bolometer signals (~30 ms), was not a limitation 

187 due to the low event rates. 

las The bolometer spectrum was divided into three energy regions: 200-400 keV, 

189 400-2000 keV, and 2000-4000 keV, as shown in Figure [6] The background rate varies 

190 by several orders of magnitude over the complete spectrum; therefore, it is useful to 

191 treat the high energy region, which contains the Q-value for 0uj3j3 decay (2527.518 ± 

192 0.013 keV), separately from the lower energy regions where the background is much 

193 higher. In addition to the Ov/3/3 Q-value, the high energy region contains the 208 T1 

194 7 line at 2614.5 keV, the 190 Pt a line at 3249 keV (including nuclear recoil), and 

195 an approximately constant background from 3-4 MeV, which is believed to be due 
we to degraded alphas. This region may also have a cosmic ray component, and was 

197 therefore investigated with this measurement. 

198 In the limit of low rates, the rate of "accidental" coincidences between muon 

199 events and bolometer events is given by: 

-^accidental = 2 ■ -Rbolo ' Rfj, ' AT (3) 

200 where -Rbolo is the bolometer event rate, i? M =4.01 mHz is the muon rate, and AT=50 

201 ms is the width of the coincidence window. Multiplying this rate by the total live 

202 time gives the expected number of accidental coincidences, which is compared to the 



number of measured coincidences in Figure 7(a) This figure shows a statistically 
significant correlation between events in the muon detector and the bolometers. 

The usual CUORICINO 0v(3(3 analysis includes an anti-coincidence cut which 
excludes any bolometer event that occurs within 100 ms of any other bolometer 
event. The bolometer anti-coincidence condition is used to reduce background, since 
the 0u/3/3 signal is expected to appear only in one bolometer. Limiting the analysis to 
single-bolometer events, the number of coincidences between the muon and bolometer 
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500 1000 1500 2000 2500 3000 3500 4000 

Energy (keV) 



Figure 6: Energy spectrum of the CUORICINO background showing the division of energy regions 
used in the analysis. No bolometer anti-coincidence cut has been applied. 



210 events is consistent with the number of expected accidentals, as shown in Figure 7(b) 

2n Evidently, the bolometer anti-coincidence cut is very effective at eliminating potential 

212 muon-induced backgrounds. 

213 The numbers of expected accidental and measured coincidences shown in Fig- 

214 ure 7(b)| provide an upper limit on the muon-induced contribution to the CUORI- 

215 C1NO background. These results are summarized in Table [2j The limits were 

216 computed by using the Feldman-Cousins method [H] to obtain an upper limit, 

217 u up , on the expected number of muon-correlated signal events. This number was 

218 converted into an upper limit on the background rate, R np , in the usual units of 

219 counts/ (keV-kg-yr) as follows: 

111 

R np = z/ up • — (4) 

Jobs X AE 

220 Here, X = 3.99 kg-yr is the total exposure (active bolometer mass times live time) 

221 and AE is the size of the energy window. The error on the energy window is taken 

222 to be on the order of the energy resolution, 7-9 keV on average. The factor / D b s = 

223 13.6 ± 1.6% is the fraction of the muons producing signal in bolometers that are 

224 also observed in the scintillators. It is obtained from the simulation described in 

225 Section [4] by taking the number of muon events which hit the bolometers and the 

226 scintillators divided by the total number of generated muon events which hit the 
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Low energy events (200-400 keV) 



Medium energy events (400-2000 keV) 



High energy events (2000-4000 keV) 
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Accidental Measured 



Accidental Measured 



Accidental Measured 



(a) Expected number of accidentals vs. measured coincidences 



Low energy events (200-400 keV) 



Medium energy events (400-2000 keV) 



High energy events (2000-4000 keV) 
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1.2 






1 

0.8 
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Accidental Measured 



Accidental Measured 



(b) Accidentals vs. measured coincidences with bolometer anti-coincidence cut 



Figure 7: Comparison of the expected number of accidental coincidences and the number of observed 
coincidences between the muon detector and bolometer signals. The error bar on the accidentals 



column represents Poisson fluctuations. Figure 7(a) includes all bolometer events in the given 
energy region, while Figure |7(b)| only includes bolometer events which pass an anti-coincidence cut 
(i.e. they do not occur within 100 ms of any other bolometer event). 
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Energy 


(A) 


M 


Upper Limits (95% CL) 








z/ u p 


R u p/10-' 3 


Low (200-400 keV) 


15.3 


15 


9.0 


83 


Mid (400-2000 keV) 


16.7 


17 


10. 


12 


High (2-4 MeV) 


1.01 





2.3 


2.1 



Tabic 2: Upper Limits (95% CL) on the contribution of muon-induced events to the CUORICINO 
background. Limits were computed using the Feldman-Cousins method. (A) is the expectation 
value of the number of accidental coincidences. M is the number of measured coincidences. v up is 
a limit on the mean number of observed muon-correlated signal events, while i? up gives an upper 
limit on the rate in counts/ (keV-kg-yr). 



227 bolometers. The uncertainty in / Q b s is the dominant systematic uncertainty in the 

228 conversion from z/ up to R up ; however, this uncertainty is much smaller than the 

229 statistical uncertainty, and has therefore been neglected in computing upper limits. 

230 After applying the bolometer anti-coincidence cut, the upper limit on the muon- 

231 induced contribution to the CUORICINO background in the 0^/3/3 region of interest 

232 is 0.0021 counts/(keV-kg-yr) at 95% confidence level. 

233 In principle, a muon (or spallation neutron) could produce long-lived (Ty 2 > 50 

234 ms) radioactive isotopes which could then decay producing a delayed coincidence 

235 signal. Based on the small number of muon events and large background, we do not 

236 expect to be sensitive to this effect. Consistent with this expectation, we find no 

237 evidence of a delayed coincidence signal. However, due to the poor sensitivity and 

238 large number of potential products (each with a different half-life and decay energy), 

239 we do not set an upper limit for delayed coincidences with the present data. 

240 5.1. Simulation Results 

241 The analysis of the simulated events was carried out in the same way as for 

242 the actual measurements. The spectrum of muon events in the various scintillators 

243 appears to be correctly reproduced in simulations. The spectrum of bolometer events 

244 was divided into the same three energy regions: 200-400 keV, 400-2000 keV, and 

245 2000-4000 keV. 

246 In Table [3j the simulated rates of bolometer events in coincidence with the muon 

247 detector are reported and compared with data (with and without imposing a bolome- 

248 ter anti-coincidence cut). The measured rates are reported after the subtraction of 

249 the expected background from accidental coincidences. In Table [4], the simulation 

250 results are reported for the total muon-induced background rate in CUORICINO. In 
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Simulation 


Measurement 




10~ 3 (counts/ (keV-kg-yr)) 


10~ 3 (counts/ (keV-kg-yr)) 


All invents 






Low (200-400 keV) 


25.0 ±0.7 


10 ±7 


Mid (400-2000 keV) 


7.91 ±0.14 


4.2 ± 1.1 


High (2-4 MeV) 


1.71 ±0.12 


1.2 ±0.4 


With Bolometer Anti- 


-coincidence Cut 




Low (200-400 keV) 


1.84 ±0.19 


< 11 


Mid (400-2000 keV) 


0.66 ±0.04 


< 1.6 


High (2-4 MeV) 


0.08 ±0.03 


< 0.29 



Table 3: Simulated and measured rates of bolometer events in coincidence with the muon detector. 
Only statistical errors are quoted. Systematic uncertainties are discussed in the text (Sections [4] 
and [5). 



Energy Total Anti-coincidence 

10~ 3 (counts/ (keV-kg-yr)) 10~ 3 (counts/ (keV-kg-yr)) 

Low (200-400 keV) 184.9 ± 1.9 7.9 ± 0.4 

Mid (400-2000 keV) 58.1 ± 0.4 3.58 ± 0.09 

High (2-4 MeV) 12.6 ± 0.3 0.53 ± 0.06 

Table 4: Simulated contribution of muon-induced events to the CUORICINO background. Only 
statistical errors are quoted. Systematic uncertainties are discussed in the text (Section [4]). 
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251 the energy region immediately surrounding the Q- value (2507.5-2547.5 keV), 

252 a value of (17.4 ± 1.3) xlCT 3 counts/(keV-kg-yr) was obtained for background in- 

253 duced by muons without any anti-coincidence cut applied and a value of (0.61 ± 

254 0.25) xl0~ 3 counts/(keV-kg-yr) with the bolometer anti-coincidence cut. 

255 6. Conclusions 

256 The bolometer anti-coincidence cut in CUORICINO appears to be a very effective 

257 tool for eliminating muon-induced backgrounds. With this cut, the measured rate 

258 of muon-correlated, single-bolometer background events was consistent with zero, 

259 and an upper limit of 0.0021 counts/ (keV-kg-yr) (95% C.L.) in the Oz//3/3 region of 

260 interest was obtained. 

261 The results of the measurement have been compared with a detailed Geant4 

262 simulation. Although the sensitivity of the experiment was not sufficient to perform 

263 a rigorous validation, the results of the measurement and simulation were generally 

264 compatible. 

265 The rate obtained for the muon-induced contribution to the CUORICINO back- 

266 ground, by measurement or simulation, is small compared to the total CUORICINO 

267 background rate of ~0.2 counts/ (keV-kg-yr) in the region of interest. Muon interac- 

268 tions also do not appear to contribute significantly to the background rate between 

269 3-4 MeV. 

270 The muon-induced backgrounds may not scale directly from CUORICINO to 

271 CUORE because of differences in the detector and shield geometry, materials, and 

272 anti-coincidence efficiency. For that reason, a detailed simulation, similar to that 

273 described in Section |4| has also been performed for muons and other external back- 

274 grounds in CUORE [19]. However, omitting subtle changes, the muon-induced back- 

275 ground rates in CUORE should be of a similar order of magnitude to those obtained 

276 for CUORICINO. The CUORE goal for the total background rate in the region of 

277 interest is 0.01 counts/ (keV-kg-yr), and both the measured and simulated values for 

278 the muon-induced background in CUORICINO are well below the CUORE goal. 
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